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We study the domain of validity of a Schwinger-Dyson approach to non-equilibrium dynamics 
when there is explicit broken symmetry. We perform exact numerical simulations of the one- and 
two-point functions of a single component A0 4 field theory in 1+1 dimensions in the classical field 
theory case, for initial conditions where {4>{x)) 7^ 0. We compare these results to two self-consistent 
truncations of the SD equations which ignore three-point vertex function corrections. The first 
approximation, which we called the bare vertex approximation, sets the three-point function to one. 
It gives a promising description for (4>(x)) = <f>(t). The second approximation, which ignores higher 
in 1/N corrections to the 2PI generating functional is not as accurate for (j)(t) for the case N = 1. 
Both approximations have serious deficiencies in describing the two-point function when cj>(Q) > 0.4. 

PACS numbers: 11.15.Pg,03.65.-w,11.30.Qc,25.75.-q 



I. INTRODUCTION 

Recently there has been much effort in finding approx- 
imation schemes to study the dynamics of phase transi- 
tions that go beyond a mean field theory (leading order in 
large- N) approach. This is an important endeavor if one 
wants a first principles understanding of the dynamics 
of quantum phase transitions. In a previous set of pa- 
pers, we have introduced a composite field x an d studied 
the validity of a next-to-leading order (NLO) truncation 
scheme of the Schwinger-Dyson (SD) hierarchy of equa- 
tions, in quantum mechanics Jl]] as well as in 1+1 dimen- 
sional field theory 0- We called this energy-conserving 
approximation, the bare vertex approximation (BVA). A 
parallel set of investigations by Berges et. al. || [|, || ff), 
have looked at a related NLO approximation based on the 
two-particle irreducible formalism (2PI-1/N), which is 
also energy-conserving. These investigators have pointed 
out that when there is explicitly broken symmetry, the 
BVA contains terms not included in the 1/N resumma- 
tion at next-to- leading order j^j. 

In this rapid communication, we present the first nu- 
merical simulations of the dynamics for a single field 1+1 
dimensional </> 4 field theory for explicit symmetry break- 
ing, where a new scale is set by initial values of the field. 
We compare the BVA with the 2PI-1/N for the classical 
field theory case, where exact comparison with numeri- 
cal simulations are possible. The differences between the 
classical and quantum field theory, are discussed in detail 
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in Refs. [^|, @], and the unbroken symmetry version of this 
paper is described in detail in Ref. Jj|. Here, we focus on 
the difference between these two approximations when 
compared with exact simulations. We also present a new 
algorithm for solving the SD equations for the broken- 
symmetry problem. 

The long-term goal of this work is finding approxima- 
tion schemes which are accurate at the small values of 
N relevant to the case of realistic quantum field theories. 
In order to maximize the possible differences between 
approximation schemes and the exact result, we choose 
to study the O(N) model for N=l. Based on our pre- 
vious studies of the quantum mechanical version of this 
model Q| we expect that by increasing N the differences 
will diminish. This is due to the fact that the SD formal- 
ism is related, but not identical, with approximations 
based on the large N expansion. As such, terms which 
may constitute next-to-next-to- leading order (NNLO) for 
the 2PI-1/N expansion are in fact NLO in the BVA. 

In our previous work g] we have shown that for the 
symmetric case, the BVA and 2PI-1/N gave results of 
similar quality. However, the agreement may have been 
misleading as in that study we only dealt with the re- 
strictive case of zero initial field expectation values. That 
choice was made solely for the purpose of reducing the 
computational cost of the calculation. In the symmetric 
case, the restrictions imposed by the choice of initial con- 
ditions, significantly simplifies the equations we need to 
solve, since many terms vanish identically by construc- 
tion. In this paper we report the first results of the full 
calculation. 

Since the 2PI-1/N expansion is a subset of the BVA, 
we can use the same code in order to generate predic- 
tions for both approximation schemes. The differences 
between the BVA and 2PI-1/N, and between these ap- 



proximations and the exact result are non-trivial. We 
will show that the BVA gives a better description of the 
order parameter <j>(t) than the 2PI-1/N. However both 
approximations suffer from deficiencies when describing 
the equilibration time of the two-point function when the 
initial value of the order parameter (f>(0) exceeds a cer- 
tain fraction of the mass parameter ji. Having established 
the efficacy of the BVA for describing the order parame- 
ter (<p(x, t)), we intend to apply this method to the 0(4) 
linear sigma model in 3+1 dimensions to study the dy- 
namics of chiral symmetry breaking in a realistic effective 
field theory. 



II. EQUATIONS OF THE 2PI EFFECTIVE 
ACTION APPROACH 



The truncation of the SD equations we are proposing 
can be obtained from a 2PI effective action P, [To| p"l|. 
Other approaches leading to these equations are found 
in (HQ. Using the extended fields notation, (f> a (x) — 
[x(xj, <f>i{x),(}>2(x),. 



(x)] , the effective action can 



be written as: 



T[<jy a , G] = S cl [<t> a ] + ^TrlnG- 1 + '-TvG^G + T 2 [G] , 

where T 2 [G] is the generating functional of the 2-PI 
graphs, and the classical action 



5',: 



i[^oi] = J dx{ ~^>i(x) [□ + x(x)] <t>i{x) 

'-x( x )}, (1) 



x 2 0) _ ^_ 

2ff ~ 9 



is that of A</> 4 field theory written in terms of the auxiliary 
field x appropriate for studying 1 /N expansions Q . The 
integrals and delta functions Sc(x,x') are defined on the 
closed time path (CTP) contour, which incorporates the 
initial value boundary condition Jl3|, |TJ, [l5| Here 
g = X/N ', and A is held fixed if we are studying the large- 
N limit. The approximations we are studying include 
only the two- loop contributions to T 2 - 

The Green function G^o[x](x, x') is defined as follows: 



S 2 S 



5cj) a {x) S(j)/s(x') 

D^(x,x') KoHx,x') 
K^{x,x') GzUx,3>) 



(2) 



where 



1 



D Q 1 (x 1 x') = — S c (x,x') , 
9 

Golj [x]( x , x') = [□ + x(x)] S^Sdx, x') , (3) 
and the off-diagonal elements are 

KqI[4>]{x,x') = Kat[4>]{x,x') = ^{x)5 c {x,x') . 



The exact Green function G a p[j](x, x') is defined by: 

S 2 W[j] 



G a p\j](x,x') = 



(4) 



8j a (x) 5jp[x') 

D(x,x') Kj(x,x') 
Ki(x,x') Gij(x,x r ) 



The exact equations following from the effective action 
are: 



[n + x (x)]cf> i (x)+K i (x,x)/i = 0, (5) 
9 
2 



X(x) = fJ- 2 + 9 - [<ft{x) + G u (x, x)/i] 



and 



where 



G a [)( x i x ') = G olp( x ' x ') + ^ais(x,x') , (6) 



, is _ ( n(x,x') CLj(x,x') 
2 ST 2 [G] 



i 6G a f3{x 1 x') 



(7) 



In the BVA, we keep in ^[G] only the graphs shown in 
Fig. [I], which is explicitly 



r 2 [G] = -- J J dxdy [G lJ {x,y)G jl {y,x)D(x,y) 

+ 2K i {x,y)K j {x,y)G ij (x ) y)] . (8) 
The self-energy, given in Eq. (Q), then reduces to: 



n(x, x ) — ^ G ran (x, x ) G mn (x, x ) , 



(9) 



(X)X ) i K rn (x, x ) G mz " (x, x ) , 
(x, x ) — i G 2m (x, x ) K m (x, x ) , 
Ey (x, x 1 ) — i [Gij (x, x') D(x, x') 4- Ki(x, x') Kj(x, x')l . 

Throughout this paper, we use the Einstein summation 

G K 

G 




G T K 

FIG. 1: Graphs included in the 2PI effective action F2- 

convention for repeated indices. 

As discussed in detail in Ref. @, the second graph in 
Fig. |l| is proportional to 1 /N 2 and is ignored in the 2PI- 
1/N approximation. Since this term is not present when 
symmetry is preserved, a comparison of these two ap- 
proaches requires studying the symmetry breaking case, 
which is done for the first time here. 
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III. UPDATE EQUATIONS FOR THE GREEN 
FUNCTIONS 



We notice from the definitions of the matrices repre- 
senting G a p(x,x') and G~l(x, x'), that the matrix ele- 
ments are not inverses of one another, but instead satisfy 
schematically: 



D- 1 D + K^ 1 K k = S e , 
D- 1 K j + K^ 1 G kj = Q, 
K i lD + G a }K k = , 
K^Kj + c;,, ; C/,v, = SijSe 



(10) 



Inverting Eq. (||), we find: 

D{x,x') — —g5c(x,x') 



(11) 



i J dxi U'(x, Xi) D(x±, x') , 



Gij(x,x') = G 0i j(x,x')Sij - / dx 1 G ik{x,x 1 ) (12) 



and 

Ki(x, x') 



X j dX2"E' kl (xi,X2)Glj(X2,x') 



J^d Xl [Do 1 +n] \x,xx) 



(13) 



x J dx 2 [K a l +Q, k \(x 1 ,x 2 )Gki(x2,x') , 
with the modified polarization 
U'{x,x') = U{x,x')- J^dx x [k^+€i k ]{x,x{) (14) 
x J dx 2 [Gqm + Sfei] ~\xi, x 2 ) [KqI + Qi] (x 2 ,x) , 
and self-energies 
X' ik (x,x') =Z ik (x,x') - [K~l +Q i ](x,x 1 ) (15) 

x J dx 2 [Do 1 + n] ~\x u x 2 ) [K- 1 + Sl k ] {x 2 ,x') . 



These update equations must be solved in conjunction 
with the one-point functions, Eqs. (||). In the approach 
of Ref . [0 , the same equations apply with the restriction 
that fli = and in the equation for Y*' ik in Eqs. ([) 15), 
the KiK k term is absent. 



IV. NUMERICAL RESULTS AND 
CONCLUSIONS 

We ran simulations for a single component (N=l) 
A (f) A classical field theory, comparing our results to ex- 
act Monte Carlo classical calculations, as described in 




FIG. 2: Results for the symmetric case, <j)(0) = 0(0) = 0. 



Rcf . [[i~7| . The classical equations become parameter free 
with a rescaling of variables as discussed in Ref. |Q. In 
particular x, t — * /ix,fit and 4> — > \/3g / '^ 2 4>- Thus by 
choosing g = 1/3, then <p and time are measured in units 
of /i and For our initial conditions, we took a ther- 
mal gaussian density matrix centered about 4>{x) = 0o 
and 4>(x) = 7To, thus we are only considering spatially 
homogeneous evolutions of <f). The details of how to ini- 
tialize the Green function equations, as well as the so- 
lutions for the symmetric case, are discussed in Ref. pj. 
The numerical procedure for solving the BVA equations 
is described in Refs. |l9|, Our update procedure con- 
served energy to better than five significant figures. 

The large scale calculations presented in this paper, 
where the full set of BVA equations is numerically solved, 
are extremely cost intensive on a computer. Hardware 
constraints impose an upper limit to the extent to which 
we can follow the time evolution of the various expecta- 
tion values. In the following, we choose to describe the 
time evolution of the system until we observe the ther- 
malization of (4> 2 (t)) in the case of the exact lattice cal- 
culation. Since in the quantum case we expect thermal- 
ization to occur faster than in the classical field theory 
case, we believe that this criteria sets a useful time scale 
for future calculations in the case of realistic quantum 
field theory models. In other words, this is the minimum 
time interval for which we must be able to follow the 
time evolution, in order to make useful predictions in the 
quantum case 

In the symmetric case, the two approximations are very 
similar and lead to the results shown in Fig. |[ We notice 
that the time-evolution of the two-point function is well 
described by the BVA. 

For the non-symmetric case, we find that the time evo- 
lution of the order parameter <fr(t) is better described by 
the BVA. For the time duration required by the exact 
calculation for the two-point function to thermalize, the 
quality of the agreement is similar for all initial values of 
4>(t). Typical results are shown in Fig. |^. Here we plot 
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FIG. 3: Plot of <^>(i)> for <f>(0) = 1, and 0(0) = 1. 



the time evolution of (4>(i)) for a case when the initial 
symmetry breaking is particularly big. Despite the fact 
that both the BVA and 2PI-1/N expansion fail to ac- 
curately describe (</> 2 (i)}, the BVA reproduces the time 
evolution of (4>(t)) quite well for almost the entire interval 
of interest. 

For the two-point correlation function the situation is 
more complicated and depends on the size of the initial 
symmetry breaking order parameter. The 2PI-1/N ex- 
pansion thermalizes for all sets of initial conditions we 
studied, however the time scale for the equilibration pro- 
cess of the two-point function, is not reproduced at large 
initial symmetry breaking values (j>(0) > 0.4. The BVA 
gives an equally good agreement as the 2PI-1/N expan- 
sion at small values of 0(0). However at larger values of 
0(0), the BVA exhibits a non-exponentially decaying os- 
cillatory behavior at intermediate times, before thermal- 
ization occurs. This behavior is illustrated in Figs. |^-|6[ 
where we gradually increase the initial value of (f)(0). For 
simplicity, we take 0(0) — for these cases. We see 




FIG. 4: Plot of (0 2 (t)} for 0(0) = 0.2, and 0(0) = 0. 
in these figures that the 2PI-1/N result is qualitatively 
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FIG. 5: Plot of <0 2 (i)) for 0(0) = 0.5, and 0(0) = 0. 
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FIG. 6: Plot of (0 2 (t)) for 0(0) = 1, and 0(0) = 0. 



better than the BVA. However, it is important to note 
that for large values of 0(0), neither approximation gives 
the correct damping behavior of the two-point function. 
This is due to the fact that three-point vertex corrections 
become important at intermediate time scales. Thus we 
conclude that although the BVA determines quite well 
the time evolution of one-point function (the order pa- 
rameter here), the time evolution of the two-point func- 
tion requires a better treatment of the three-point vertex. 

The accuracy of the BVA depends on the size of the 
initial symmetry symmetry breaking. This seems to indi- 
cate that for a regime where thermalization occurs rela- 
tively slowly, such as in the case of this classical field the- 
ory model, higher order corrections do eventually become 
increasingly important at later times. The quantum field 
theory case is qualitatively different. In the quantum 
case, thermalization is expected to settle in much faster. 
Therefore, one can hope that in the quantum case the 
noise will not have enough time to buildup and conse- 
quently will result in a lesser distortion of the BVA pre- 
dictions. Future studies must go beyond the BVA and 
investigate the magnitude of these further corrections, 
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especially in the quantum case. 

In the above approximations, we set the vertex func- 
tion to its bare value of unity. However the formalism 
tells us how to determine the corrections contained in T2. 
The three-point vertex function is found from the self- 
energy S Q/3 by: 

r afiji.X 7 X , X ) 

= fapySix, x')S(x, x") + - - , (16) 

where / ,.y = /j,oj = fi,j,o = <%■ The integral equa- 
tion generated from this equation for T a 0j(x,x',x") is 
discussed in Ref. Thus one way of including vertex 
corrections is to use the Green functions calculated here 
to obtain an approximate nonlocal vertex function from 
Eq. (16), reintroduce the vertex function into the exact 
SD equation for the two-point function equation and it- 
erate until a stable solution is found. A second approach 
is to include three loop terms in T2 and solve the result- 
ing set of equations. These approaches will be explored 
elsewhere. 

To conclude, in our quest for an accurate assessment 
of NLO versus mean-field effects in quantum field the- 
ory, the present work represents an important stepping 
stone. In this paper we report the first results of the full 
BVA calculation, and in doing so we conclude work done 
earlier in a more restrictive regime Q . Our goal remains 
the study of the quantum field theory case, but the clas- 
sical field theory case is still interesting since this is the 
high temperature limit of the quantum case. We want 



to make sure that we know to what extent this limit is 
reproduced. This is even more important since this is 
the only case for which exact calculations are possible, 
and thus the classical field theory 0(1) model with one 
spatial dimension represents a key benchmark. 

Having established the domain of applicability of the 
BVA in this extreme limit, it is clear that two avenues 
must be explored in the near future: first we will apply 
the BVA formalism to the study of the quantum linear 
sigma model, and other realistic models, in 3+1 dimen- 
sions, and in doing so we will obtain first insights into the 
reliability of mean-field approaches in studies of the real- 
time evolution of realistic quantum field theory models. 
Next we will investigate the role of further corrections to 
the BVA, which will be the real test of the accuracy of 
the BVA formalism. This project will involve the devel- 
opment of practical numerical algorithms and, possibly, 
require the next generation of supercomputers. 
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